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ABSTRACT 

We analyze the optical-near infrared spectra of 33 quasars with redshifts 3.9 < z < 6.4 to investigate the properties of dust extinction 
at these cosmic epochs. The SMC extinction curve has been shown to reproduce the dust reddening of most quasars at z < 2.2; 
we investigate whether this curve also provides a good description of dust extinction at higher redshifts. We fit the observed spectra 
with synthetic absorbed quasar templates obtained by varying the intrinsic slope (a A ), the absolute extinction (A i000 ), and by using 
a grid of empirical and theoretical extinction curves. We find that seven quasars in our sample are affected by substantial extinction 
(A 30 oo ^ 0.8), and characterized by very steep intrinsic slopes (a A < -2.3). All of the individual quasars require extinction curve 
deviating from that of the SMC, with a tendency to flatten at A < 2000 A (in the rest frame of the source). However, due to the 
uncertainties in the individual extinction curves the SMC is still (marginally) consistent with the data in most cases. We obtain a 
mean extinction curve at z>4, both by performing a simultaneous fit of all quasars and by averaging the extinction curves inferred 
for individual quasars. In the case of broad absorption line quasars (which are generally more absorbed by dust and possibly in a 
younger evolutionary stage), the mean extinction curve deviates from the SMC at a confidence level > 95%. The difference between 
S_H ■ extinction curves in quasars at z>4 and in quasars at lower redshift is indicative of either a different dust production mechanism at high 

redshift, or a different mechanism for processing dust into the ISM. We suggest that the same transitions may also apply to normal, 
star-forming galaxies at z>4. In particular, the observed change in the average spectral slope of galaxies at z>4 may be partially 
ascribed to a variation in the extinction curve, rather than a lower dust content at high redshift. In this scenario, the extinction curve 
inferred at z>4 would imply a cosmic star-formation corrected for dust attenuation a factor of ~ 2 higher than estimated in the past. 
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^ ; 1. Introduction Extinction curves that differ from the SMC, LMC, and MW 

^J. , have been obtained in the case of local AGNs. In the case 

\q . Dust represents one of the key ingredients of the Universe by of mildly obscured AGNs, various authors (Gaskell & Benker 

O playing a crucial role in both the formation and evolution of the 2 007, Maiolino et al. 2001) have found indications of a very flat 

O stellar populations in galaxies as well as their observability. extinction curve in the UV. Richards et al. (2003) and Hopkins 

Extensive measurements of the local interstellar extinction e t al. (2004) found that mildly reddened QSOs (including broad 

; exist along hundreds of lines of sight in the Milky Way (MW) absorption line, BAL, QSOs) at z < 2.2 are characterized by 

| ■ tH , (ISavage & Ma this 1979), and along tens of lines of sight in the SMC-like extinction curves. The two results are not inconsistent 

^ ■ Large Magella nic Cloud (LMC) and in the Small Magellanic with each other: while in mildly reddened QSOs the line of sight 

H \ Cloud (SMC) dFitzpatrick 1 1 1 989t) . The wavelength-dependence passes through a low density and often outflowing medium ex- 

. r* , of the dust extinction and spectral features in extinction curves posed to the QSO radiation, in the case of obscured AGNs we 

are useful for constraining the size distribution of dust grains and are probably piercing through dense clouds where dust growth 

revealing the dust chemical composition. The extinction curves may lead to a flattening of the extinction curve (s imilarly to that 

which characterizing the local Universe have been studied by observed in Galactic clouds, Card elli et alj|1989t) . 
several autho rs (e . g.. ICardelli et alj|1989t iMathis et all 1 19771) . 

In particular, iPeil (119921) . by adopting a graphite-silicate grain In the case of star-forming galaxies, it is generally difficult 

model, finds that the MW reddening requires roughly an equal to infer the extinction curve since, in contrast to quasars, the 

amount of graphite and silicate grains, while the SMC extinc- simple dusty screen case does not apply, a more accurate de- 

tion curve is dominated by silicate grains (see also Weingartner scription being dust mixed with the emission sources (stars and 

& Draine 2001), and the LMC represents an intermediate situa- ionized gas). As a consequence, the observed spectrum heavily 

tion. depends on th e rela tive distribution of emitting sources and dust. 



Calz etti et alJ (1 19941) inferred an "attenuation curve" for star- 
* Based on data obtained at the VLT, through the ESO programs forming galaxies that determines the effective reddening of the 
074.B-0185 and at the Telescopio Nazionale Galileo. global galaxy spectrum relative to the case of no extinction, and 
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is flatter than empirical extinction curves (SMC, LMC, MW). 
The Calzetti attenuation law can be interpreted in terms of dust 
mixed with the emit ting sources, and with an SM C-like intrinsic 
extinction curve dInouel2005tlGordon et alJ [l997). although the 
Calzetti law may also be recovered by assuming MW-type dust 
and different geometries for the emitting sources and dust media 
dPanuzzo et al.ll2007blPierini et alJl2004 . 

At high redshift, the spectral energy distribution (SED) of 
lyman break galaxies (LBGs) is often modeled with the Calzetti 
attenuation curve, although the SMC extinction cur ve provides a 
better description for young star-forming galaxies ( Redd v~et al.l 
120101) . indicating that a screen-like dust distribution exists in the 
latter sour ces (probably d ust associated with the starburst wind). 
However, iNollet all (2009) found that more massive sources at 
z~2 display an absorption dip at 2175 A, typical of MW-like 
extinction curve, suggesting that these objects are in a more 
evolved evolutionary stage with dust properties similar to local 
disk galaxies. 

At the highest redshifts probed so far (z ~ 6), the existence 
of large amounts of dust has been inferred from mm and submm 
observations of quasars host galaxies. The far-infrared (FIR) lu- 
minosities of these objects are consistent with thermal emission 
from warm dust (T < 100 K) with dust masses M dust > 10 8 M , 
assuming that the domin ant dust heating mechanism is radiation 
from young stars (e.g. ICarilli et alJ l2Q0lb lOmont et ail [200 1 ; 
iBertoldi et ai]|2003l: iBeelen et al.ll2006HWang et al.ll2008l) . 

Many efforts have been made to explain the enormous 
amount of dust observed within 1 Gyr of the Big Bang. 

Dust formation occurs in environments characterized by 
temperatures of about 1000 - 2000 K and high density (n ~ 
10 8 - 10 15 cm~ 3 ). The proper combination of these conditions is 
mainly found in the atmospheres of low and intermediate-mass 
(M < 8 M Q ) evolved stars during the asymptotic giant branch 
(AGB) phase, from where dust is transported into the ISM by 
means of stellar winds. However, it is generally assumed that the 
timescales required for the bulk of the AGB population to evolve 
(a few 10 8 to a few 10 9 yr) are comparable with the age of the 
Universe at z ~ 6 (~ 1 Gyr); as a consequence, the possibility 
that evolved stars can explain the copious amount of interstellar 
dust observed at these redshifts has been largely questioned (e.g. 
Morgan & Edmunds 2003, but see also Valiante et al. 2009). 

Core-collapse type II supernovae (SNell) may constitute a 
valid alternative to the dust production at high redshift. They 
provide a rapid mechanism of dust enrichment since the life- 
times of their progenitors are short (~ 10 6 yr), and dust grains 
are predicted to form in their expanding ejecta a few hundred 
days after the explosion. This scenario is corroborated by the in- 
frared emission and the asymmetric evolu tion of optical lines 
observed in the spectra of nearby SNe ([W ooden et al. 1993; 
Elmhamd iet al.l2003tlSugerman et al.|2006l; lMeikle et al.[2007; 
Kotak et all 120061; IWesson et al.l l2010). as well as by infrared 



and submm observations of SN remnants dMorgan et al.l l2003: 
iRho et al.ll2008l 120091; iGomez et al.ll2009l) . These observations 
suggest that the total mass of dust produced by SNe is of the or- 
der of a few times 10~ 2 M per SN, although these estimates are 
disputed. For example, dust yields as high as 0.1 - 1.0 M /SN 
have been claimed for the SN remnants Cassiopeia A and K epler 
dMorgan et al.ll2003b IGomez et al.ll2009t iDunne et al.ll2009t) . 

From a theoretical point of view, several authors have 
proposed that the ejecta of SNell (and also pair instabil- 
ity SNe, PISNe) are an efficient dust formation mecha- 
nism at high redshifts (|Todini & Ferraral 1200 it ISchneider et al.l 
l2004t iHirashita etall 120081) . These models predict both the 
SNell/PISNe extinction curves and the masses of dust produced 



by each class of SN. For instance, SNell from progenitors with 
12-40 M are predicted to produce dust masses of about 
0.1 - 2 M /SN. The discrepancy between the mass of dust in- 
ferred from most of the observations (< 5 x 10~ 2 M ) and the- 
oretical predictions (> 0.1 M Q ) can be reconciled by includ- 
ing the effects of the reverse shock, which can destroy a sub- 
stantial (> 80%) fraction of the dust formed in the SN ejecta 
dBianchi & Schneiderll2007l:lNozawa et al.ll2007HHirashita et al.l 
120081) . 

The scenario of SN-like dust formation in the early Universe 
has been investigated observationally by studying the extinction 
curve in z~6 sources, and in particular by analyzing the red - 
dened quasar SDSSJ1048+46 at z = 6.2 (iMaiolino et al Il20 04a). 
the GRB050904 aftergl ow at z = 6.3 (IStratta et alJ2007l) and the 
GRB071025 at z ~ 5 dPerlev etal.ll2010l) . These studies indi- 
cate that the inferred dust extinction curves differ from those ob- 
ser ved at low redshif t, while being in very good agreement with 
the ITodini & Ferraral (1200 ll) predictions of dust formed i n SNe 
ejecta. However, it must be mentioned that iLi & Liangl ([2008) 
also analyzed the extinction of the GRB050904 host galaxy, 
finding a somewhat different result from that of St ratta et al.l 
(2007), but by assuming for the GRB intrinsic optical contin- 
uum a much shallower power-law slope /3 op , than that inferred 
from the X -r ay dat a (more specifically p opt ~ 0.25 vs fix ~ 1.2). 
Zafar et alJ (l2010l) revised the photometric measurements of 
GRB050904 reported in the literature and investigated whether 
this GRB is significantly absorbed by dust. 

The question of whether the SN-dust formation is efficient 
enough to account for Mdust > 10 8 M at z ~ 6 remains a topic 
of debate. lDwek et all d2007) applied the results of an analytical 
model for the evolution of the gas, dust, and metals in high-z 
galaxies to one of the most distant quasars (SDSSJ1148) at z — 
6.4. Based on the assumption that SNe are the only available dust 
suppliers at high redshifts, these authors find that their model can 
explain the observed dust mass, provided that the IMF is top- 
heavy and the dust condensation efficiency is as high as 100%; 
they also sugg est that substantial grain g ro wth may occur in the 
ISM (see also lMichalowski et al.l2010albl) . IValiante et al.l (f2009) 
are able to reproduce the dust mass observed in SDSSJ1148 at 
z = 6.4, under much less extreme conditions, by means of a 
chemical evolutionary model that includes the role of AGB stars 
as dust suppliers at high redshifts. By comparing their model 
predictions with the mass of dust inferred from SDSSJ1 148, they 
find that AGB stars can account for 50-80% of the total mass 
of dust, while SNe may account for the remaining fraction (50- 
20%), depe nding on the star-fo rmation history of the host galaxy. 

Finally, lElvis et ail (120021) present an alternative to the nor- 
mally assumed stellar origin of cosmic dust, by proposing a 
dust condensation scenario connected to quasar outflow winds. 
According to their scenario, the conditions for dust formation 
may occur within the clouds responsible for the broad emis- 
sion lines in expansion through the quasar outflowing wind (pro- 
vided that metals, and in particular refractory elements, are al- 
ready available). Whether this mechanism can explain the huge 
amount of dust at high redshift remains unclear and depends on 
the quasar outflow rate across its evolutionary history, the quasar 
nuclear metallicity evolution, and the grain condensation effi- 
ciency within the outflowing clouds. 

To summarize, the dust production mechanism at high-z is 
still a hotly debated issue. In this work, we study the properties 
of dust at 3.9 < z < 6.4 by analyzing the optical/NIR spectra of 
a sample of ~ 30 QSOs. In particular, we investigate whether 
the dust extinction inferred from high redshift quasar spectra 
can be described in terms of the SMC extinction curve, as for 
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quasars at z < 2.2 (Richards et al. 2003; Hopkins et al. 2004), 
or it reflects variations in the dust grain size distribution and/or 
in the dust chemical composition. The outline of this paper is 
as follows. Section 2 presents the observational data. In Sect. 
3, we describe the properties of the extinction curves adopted 
in this study. Moreover, we discuss the results of the analysis 
of the quasar spectra that require substantial dust reddening and 
we quantify the deviation of the mean observed extinction curve 
from the SMC. The implications of our findings are discussed in 
Sect. 4, while, in Sect. 5, we summarize the main results of our 
study. 



2. Spectroscopic data and sample 

Most of the spectra used in this study were obtained in various 
observing runs at the Telescopio Nazionale Galileo and at the 
ESO Very Large Telescope. Many of these spectra were already 
partly presented in previous works (Maiolino et al. 2003, 2004b, 
Juarez et al. 2009). The data reduction process is presented in 
these previous papers, and here we summarize only the instru- 
mental setup. The majority of the s pectra were ob tained with the 
TNG near-IR spectrometer NICS ( Baff a et al .1200 lb by using the 
AMICI prism, which allows us to observe the entire 0.9-2.3 //m 
spectral range in a single observation, though at low spectral res- 
olution (R~50-75, which is enough for the purposes of inves- 
tigating the dust extinction through the spectral shape). A few 
quasars were also observed with NICS at intermediate disper- 
sion (R~300) through the IJ and HK grisms (mostly to more 
clearly resolve a few absorption and emission line features). 

We also include a number of quasars at z<5 observed with 
FORS2 by using the grism GRIS150I, which covers the 0.7- 
1.0 //m range for R~300. However, in most of these cases the 
limited rest-frame spectral coverage does not allow us to con- 
strain the dust reddening. 

The selection criterion for quasars at z>5 was simply to in- 
clude all quasars known at the epoch of the observations and 
visible during the observing runs. In the redshift range 3.9<z<5, 
we gave preference to known BAL quasars, since these tend to be 
more often affected by dust extinction, hence more suitable for 
investigating the extinction curve. As a consequence, the large 
number of BAL quasars in our sample of quasars at 3.9<z<6.4 
simply reflects our selection criterion. 

We also include in our sample the near-IR spectra of five 
z~6 quasars by Jiang et al. (2007), which were obtained with a 
spectral resolution R ~ 800 at GEMINI South by using GNIRS, 
which provides a simultaneous wavelength coverage from 0.9 to 
2.5 [an. Finally, we also include the NIR spectrum of CFHQS 
J 1509- 1749 (z = 6.12) obtained by Willott et al. (2007) at 
GEMINI South by using GNIRS and the NIR spectrum of ULAS 
J0203+0012 (z = 5.72) by Mortlock et al. (2009) at GEMINI 
North with GNIRS. 

All the quasars adopted in our analysis are reported in Table 
1 along with important information about the spectroscopic 
data. The spectra of most quasars are shown in electronic form 
(Appendix A), except for the reddened targets (the most inter- 
esting cases), which are shown in the next section. 

Hereafter, we refer to each object by using only the acronym 
of the survey and the first four digits of its right ascension and 
declination. 



3. Determining the extinction curve at high redshift 

3.1. The method 

The spectrum of each quasar is fitted with the following expres- 
sion 

F A = CF'\ -\ 10 25 "3000 , (1) 

\3000A/ 

where C is a normalization constant, F' A is a quasar template 
spectrum, ax the intrinsic slope of the unreddened spectrum, 
A 3000 the absolute extinction at 3000 A, and A\ /A 3000 is the ex- 
tinction curve normalized at 3000 A, which is derived from the 
library of empirical and theoretical extinction curves discussed 
in the next section. In this study, we adopt the template com- 
piled by Reichard et al. (2003), obtained by considering 892 
quasars classified as non-BAL, whose continuum is described 
by a power-law with index a^, = -1.62 (F' A oc A aAj ). Therefore, 
the term (/l/3000) (I 62+ff ,) allows us to force the slope of the tem- 
plate to the value a^. We allow ax to vary in the interval [-2.9; 
-0.2], which is the range encompassed by more than 95% of the 
QSOs in the SDSS sample of Reichard et al. (2003). The nor- 
malization of the intrinsic quasar spectrum is allowed to vary, 
with the only constraint being that the intrinsic quasar bolomet- 
ric luminosity should not exceed the Eddington value, within a 
factor o f 1.5, to allow for un certainties in the bolometric cor- 
rection (Risali ti~& Elvis Il2004l) . For each quasar, the Eddington 
luminosity is inferred from black hole masses obtained from the 
literature (Willott et al. 2003; Jiang et al. 2006; Kurk et al. 2007), 
or estimat ed from our own spectra by us ing the virial calibrations 
given in lVestergaard & Peterson! (120061) . 

In the fitting procedure, we avoid spectral regions affected 
by the following emission lines: Lyc+NV [1215.67; 1280] A; 
Ol+Sill [1292; 1320] A; SilV+OIV] [1380; 1418] A; CIV 
[1515; 1580] A; A1III+CIII] [1845; 1945] A. We also exclude 
the region 2210-3000 A, since it is generally characterized by a 
prominent Fell bump. Moreover, in the case of BAL QSOs, we 
identify, and discard from the fitting, the associated absorption 
troughs. Spectral regions with some instrumental artifacts were 
also removed. Finally, we avoid the following regions in the ob- 
served frame, which are affected by poor atmospheric transmis- 
sion: 1.33-1.5 ym\ 1.8-1.95 fan; > 2.3 /urn. After selecting the 
spectral regions to be included in the analysis, we rebin the ob- 
served spectra to a resolution R ~ 50, which is about the spectral 
resolution delivered by the Amici prism in most NICS-TNG ob- 
servations. However, we checked that the results do not depend 
on the binning adopted. 

For each quasar observed in our sample, we find the best fit- 
ting ax, A3000, and the extinction curve Ax that minimize the^- 2 . 
Errors in our spectra (especially in the near-IR) are dominated 
by calibration uncertainties. In particular, the intrinsic shape of 
the telluric calibration standard star (used also to calibrate the 
relative instrumental response) is generally known with an un- 
certainty of about err re i = 5%, which is therefore assumed as 
the average uncertainty in our spectra. In some cases, some re- 
binned spectra regions are affected by a noise higher than 5%; in 
these cases, we rebinned the spectra again to reach noise levels 
below the calibration uncertainty of 5%. 

Besides the best fit x 2 Of hf)> we compute the x 2 for the 
case where no extinction is assumed (i.e. A 3000 = 0), defined as 
X%oexf ^ e c ^ m tnat a gi ven quasar spectrum requires extinc- 
tion if the probability of getting a^- 2 lower than^ . is > 90% 
(i.e. P(Xfj oe S) ^ 90%), and, at the same time, if the probability 
of getting ax 2 lower than^ 2 F is < 68% (i.e. P(x 2 BF ) < 68%). 
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Table 1. Sample of high-z quasars 



Name 


z 


Ref. 


Instrument 


Dispersed 


SDSS J000239.4+255035 


5.80 




NICS-TNG 


Amici+IJ 


SDSS J000552.3-000656 


5.85 




NICS-TNG 


Amici 


SDSS J001714.66-100055.4 b 


5.01 


HI 


FORS2-VLT 


GRIS 1501+27 


SDSS J012004.82+141108.2 b 


4.73 


[1] 


NICS-TNG 


Amici 


SDSS J015642.11 + 141944.3 b 


4.32 


HI 


NICS-TNG 


Amici 


ULAS J020332. 38+001229. 2 b 


5.72 


[4] 


GST 


GNIRS 


SDSS J023 137.6-072855 


5.41 




FORS2-VLT 


GRIS 1501+27 


SDSS J023923.47-081005.1 b 


4.02 




NICS-TNG 


Amici 


SDSS J033829.3+002156 


5.00 


HI 


FORS2-VLT 


GRIS 1501+27 


SDSS J075618. 1+410408 


5.07 


HI 


NICS-TNG 


Amici 


APM 08279+5255 b 


3.90 




NICS-TNG 


Amici 


SDSS J083643. 8+005453 


5.80 




GST 


GNIRS 


SDSS J085210.89+535948.9 b 


4.22 


HI 


NICS-TNG 


Amici 


SDSS J095707.67+061059.5 


5.16 




FORS2-VLT 


GRIS 1501+27 


SDSSp J1021 19. 16-030937.2 


4.70 




FORS2-VLT 


GRIS 1501+27 


SDSS J103027. 1+052455 


6.28 




GST 


GNIRS 


SDSS J104433.04-012502.2 b 


5.78 


[2] 


GST 


GNIRS 


SDSS J104845.05+463718.3 b 


6.20 


HI 


NICS-TNG 


Amici+IJ +HK 


SDSS Jl 14816.6+525150 


6.40 


HI 


NICS-TNG 


Amici 


SDSS J120441.7-002150 


5.05 


HI 


FORS2-VLT 


GRIS 1501+27 


SDSSp J120823. 8+001028 


5.27 




FORS2-VLT 


GRIS 1501+27 


SDSS J130608.2+035626 


5.99 


11 1 


GST 


GNIRS 


SDSS 1141111.3+121737 


5.93 


| 1 


GST 


GNIRS 


CFHQS J150941.78-174926.8 b 


6.12 


| -}| 


GST 


GNIRS 


SDSS J 160254.2+422823 


6.07 




NICS-TNG 


Amici+IJ 


price TI/^AQOn QO i llld C 

oDss J 1oU3zU.89+0 /21U4.J 


4.39 


m 


JNrCo-lJNlj 


Amici 


SDSS J160501. 21-01 1220. 6 b 


4.92 




NICS-TNG 


Amici+IJ 


SDSS J161425. 13+464028.9 


5.31 




NICS-TNG 


Amici 


SDSS J162331. 8+311201 


6.26 




NICS-TNG 


Amici 


SDSS J162626.50+275132.4 


5.20 




NICS-TNG 


Amici 


SDSS J163033.9+401210 


6.06 




NICS-TNG 


Amici 


SDSS J220008.7+001744 


4.77 




NICS-TNG 


Amici 


SDSS J22 1644.0+00 1348 


4.99 




NICS-TNG 


Amici 



"Resolution: Amici (R ~ 50 - 75), IJ (R ~ 300), GRIS (R ~ 300), GNIRS (R ~ 800), HK (R ~ 300); b BAL quasars; [1] This work and Juarez et 

al. (2009); [2] Jiang et al. (2007); [3] Willott et al. (2007); [4] Mortlock et al. (2009) 



It should be noted that the method to infer the extinction 
curve adopted in this paper, i.e. by fitting template spectra 
and extinction curves to the observed spectra, differs from the 
method employed in Maiolino et al. (2004a), who inferred di- 
rectly the extinction curve by inverting Eq. 1 . Although the lat- 
ter method has the advantage of not assuming any shape for the 
extinction curve a priori, it does require assumptions about the 
intrinsic spectrum of the QSO, in terms of both normalization 
and slope; in particular, in Maiolino et al. (2004a) a quasar lu- 
minosity close to the Eddington limit was assumed and spectral 
slopes varying within an interval narrower than adopted here. In 
contrast, with the "fitting method" used in this paper, the intrin- 
sic spectrum of the QSO, and in particular its normalization and 
slope, are quantities that are fitted by the procedure. 

3.2. Library of extinction curves 

To investigate the evolution of the dust properties across cos- 
mic time, we consider a grid of extinction curves to charac- 
terize the extinction produced by dust in the rest-frame wave- 
length range 1000 A < A < 5000 A. We first consider the 
empirical curves that describe the dust extinction in the local 
Universe (ISavage & Mathisl 1 1 979t lFitzpatrickl [l989). i.e., the 
SMC, LMC, and MW extinction curves (Fig. Q] upper left hand 
panel). We then use the extinction curves expected by Typ e 
II SNe dust models as predicted by Todi ni & Ferraral ([2001), 



Bian chi & Schneiden(l2007l). and Hirashita et al. (2008). The ex- 
tinction curves of iTodini & Ferraral d2001l) . shown in the right 
hand upper panel of Fig. [TJ depend on the metallicity of the SNe 
progenitors; the right-hand upper panel of Fig.[T]shows the cases 
of Z = (TF1, magenta solid line), Z = Z Q (TF2, green dotted 
line), and Z = 10~ 2 Z o (TF3, red dotted-long-dashed line). In 
their model, iBianchi & Schneider! (120071) take into account the 
possibility that the freshly formed dust in the SNe envelopes 
can be destroyed and/or reprocessed by the passage of the re- 
verse shock. They therefore follow the evolution of the dust con- 
densed from the the time of formation in the ejecta to its sur- 
vival in the SNe remnants a few thousand years later. The right- 
hand upper panel of Fig. Q] shows the extinction curve predicted 
by the model before (BS1, cyan dotted line) and after (BS2, or- 
ange dotted-short-dashed line) the reverse shock, assuming solar 
metallicity for the SNe p rog enitor. Both the extinc t ion cu rves by 
ITodini & Ferraral d2001l) and lBianchi & Schneiderl (120071) are av- 
eraged over the Salpeter initial mass function. 
In the bottom l eft-hand panel ofTtg.Q] we show the results of 
the models by Hira shita et al.l d2008l) for dust produced by SNe 
of progenitor mass 20 M . In this work, the effect of the reverse 
shock is also consi dered, for various dens ities of the ambient 
medium. Moreover, iHirashita et al.l ((2008) study two extreme 
cases for the mixing of the elements constituting the dust grains: 
the first is the unmixed case, where the original onion-like struc- 
ture of elements is preserved in the helium core; the second is the 
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Fig. 1. Upper left panel: Empirical curves for dust extinction in the local Universe (Savage & Mathisl 1 1 979b lFitzpatrick1 [T989). 
Upper right panel: Theoretical extinction curves predicted by Type II SNe dust models proposed by Todini & Ferrara (2001) and 
Bianchi & Schneider (2007). Lower left panel: Results of the Type II SNe dust model by Hirasita et al. (2008). Lower right panel: 
Predictions for extinction curves produced by PISNe dust (Hirasita et al. 2008). The main properties of the theoretical models shown 
in this figure are mentioned in Sect. 13.21 



mixed case, characterized by a uniform mixing of the elements. 
The solid magenta (dashed cyan) line shows the mixed case, af- 
ter the reverse shock, for a medium characterized by hydrogen 
number density of n# = 1 crrT 3 («# = 0.1 cm -3 ). Hereafter, we 
label these models as HI and H3, respectively. The dotted green 
line represents the unmixed case for freshly formed dust grains 
(H2), while the dotted-dashed orange line shows the predictions 
for the unmixed case after the reverse shock, for = 0.1 cirT 3 
(H4). 

iHirashita et alJ d2008l) also investigated the case in which the 
extinction is caused by dust produced by pair instability super- 
novae (PISNe) with a progenitor mass of 170 M . In the right- 
hand bottom panel of Fig. [TJ we show the results of their anal- 
ysis for a medium characterized by n# = 1 cirT 3 in the mixed 
(PISN1, solid magenta line) and unmixed (PISN2, dotted green 
line) cases. 

We note that our analysis does not really aim to distinguish be- 
tween different SN-dust progenitors using the inferred extinc- 
tion curve. We more simply aim to adopt a library of extinction 
curves that are physically plausible. 



Predictions for the extinction curve produced by dust formed 
in quasar winds (Elvis et al. 2002) have not yet been developed, 
hence this scenario cannot be included in our analysis. 

Finally, we note that the dust that causes the reddening 
of a quasar spectrum may not be produced by a single ex- 
tinction curve among those discussed above, but by a mixture 
of dust grains with different properties, origins, and extinction 
curves. Therefore, besides the grid of extinction curves men- 
tioned above, we also consider mixed extinction curves. In par- 
ticular, since the dust extinction inferred from quasar spectra at 
< z < 2 can be described in terms of the SMC extinction curve 
(Richards et al. 2003; Hopkins et al. 2004), we only mix the 
modeled extinction curves with the SMC, which is our low-z ref- 
erence for quasars. For exam ple, in the case of zero- metallicity 
extinction curve predicted by Todini & Ferrara ( 200 1 ) for a SNe- 
like dust, we consider the following family of extinction curves 

At = {\-p)A s A MC + P A T /\ (2) 

where < p < 1 defines the fraction of each SNe-type con- 
tributing to the extinction curve. We consider variations in the p 
parameter with steps of Ap — 0.1. 



5 



S. Gallerani et al.: The extinction law at high redshift and its implications 
Table 2. Best-fit model parameters for reddened quasars 



Name 


z 




"A.BF 


^3000,Sf 


&A,Noext 


a A,SMC 


P(X BF ) 


P^Noext^ 




d.o.f. 


SDSS J1148 


6.40 


70% SMC+30%PISN2 


-2.90 


0.82 


-1.75 


-2.57 


0.43 


0.99 


0.53 


14 


SDSS J1623 


6.26 


60% SMC +40% PISN2 


-2.63 


0.92 


-1.39 


-2.04 


0.44 


0.94 


0.47 


12 


SDSS J1048 a 


6.20 


30%SMC + 70% H4 


-2.90 


1.43 


-1.75 


-2.90 


0.14 


1.00 


1.00 


15 


ULAS 0202 a 


5.72 


50% SMC+50% PISN2 


-2.55 


1.69 


-0.73 


-1.46 


0.40 


0.99 


0.63 


16 


SDSS J1626 


5.20 


60% SMC+40% PISN2 


-2.90 


1.64 


-1.00 


-2.36 


0.08 


1.00 


0.23 


24 


SDSS J0852 a 


4.22 


40% SMC + 60% PISN2 


-2.33 


2.00 


-0.50 


-1.14 


0.01 


1.00 


0.92 


26 


APM0827 3 


3.90 


40% SMC+60% PISN2 


-2.67 


1.94 


-0.90 


-1.55 


0.00 


1.00 


0.04 


30 



11 BAL quasars. 

Table 3. Simultaneous best-fit model parameters for reddened quasars 



Sample 


AlBF 




P^Noext^ 


P<Xiuc) 


d.o.f. 


ALL 


60% SMC+40%PISN2 


0.10 


1.00 


0.61 


161 


BAL 


40% SMC+60% PISN2 


0.17 


1.00 


0.94 


99 


noBAL 


70%SMC + 30% PISN2 


0.02 


1.00 


0.08 


58 



As far as this work is concerned, we take into account empir- 
ical and theoretical extinction curves, while we do not consider 
any kind of parametric curve (e.g. Li & Liang 2008). We found 
that the use of parametric curves is hampered by a strong degen- 
eracy among the various free parameters, normalization constant 
C, the continuum slope a% and the absolute extinction A3000, 
which enter in Eq. 1 . We note that this problem is alleviated in 
the case of the GRB afterglow spectra, for which the intrinsic 
slope of the continuum can be inferred from X-ray observations, 
which are not affected by dust. 

3.3. Extinction curves of individual quasars 

We find that seven quasars require significant dust extinction. 
We note that the large fraction of extincted quasars found by 
ourselves (~ 23%) is much higher than the one found in lower 
redshift samples (~ 5%, Hopkins et at. 2004, Reichards et al. 
2003). However, this fraction does not reflect the true fraction of 
reddened quasars at high redshift, being simply a selection ef- 
fect, because we preferentially selected BAL quasars, which are 
generally affected by higher extinction than non-BAL quasars. 

In Table|2] we report the quasars that according to our analy- 
sis are reddened, along with the extinction curve minimizing the 
X 2 , described in Sect. 13.21 the best-fit model parameters ctx,BF 
and A3ooo,bf, the slope in the "no extinction" case a A ^ oext , the 
probability associated with the best-fit values P(x 2 BF ), the "no 

extinction" case P(Xff oext ), the "pure SMC" case P(Xsmc^ anc * 
the degrees of freedom (d.o.f.) associated with each spectrum. 

In Figs. [2]|3] we show the observed/rebinned spectra (black 
solid line/blue filled squares), and the best-fit spectra (red dashed 
line) for these reddened quasars. Figure [4] (left panel) shows 
the resulting best-fit extinction curves for each reddened spec- 
trum (BAL quasars with solid lines and non-BAL quasars with a 
dashed lines). 

The inferred extinction curves for all these reddened quasars 
at z>4 deviate from the SMC, with the tendency to flatten at short 
wavelengths (A RF < 2000 A) relative to the SMC. However, 
the SMC extinction curve, although not providing a best fit for 
any object, is still marginally consistent with the data at the lcr 
level in most of the reddened quasars. Only for SDSSJ1048 and 
SDSSJ0852 is the SMC extinction curve ruled out with high con- 
fidence. 



Finally, we note that three of the seven reddened quasars re- 
quire an extreme slope of a = -2.9. While in Sect. I3.5l we dis- 
cuss selection effect that probably favor the selection of high-z 
reddened quasars with very blue slopes, it is worth checking the 
effect of restricting the allowed spectral slopes to a more "rea- 
sonable" range, more specifically a > -2.6, which still charac- 
terize 90% of quasars at lower redshifts. The result of changing 
the maximum slope is shown in Fig. |4] (right hand panel). Only 
for SDSSJ1048 does the shape of the extinction curve change, by 
preferring TF1 (as in Maiolino et al. 2004a), but the differences 
from previous cases are not large. We also restricted the spectral 
slope range to [-2.0; -1.2], which is the interval encompassed 
by 68% of quasars at lower redshifts, finding that the average 
extinction curve (see next section) does not change significantly. 

3.4. The average extinction curve at high redshift and the 
deviation from the SMC 

To determine how significantly the extinction curve deviates 
from that of the SMC at z>3.9, and provide an average descrip- 
tion of the extinction curve at high redshift, we follow two ap- 
proaches. 

We first compute the mean of the inferred extinction curves 
to provide an empirical extinction law at z > 3.9. The result is 
shown in Fig. [3 where the black solid line represents the mean 
of the inferred extinction curves. Hereafter, we refer to this 
extinction curve as the "mean extinction curve" (MEC). The 
shaded area illustrates the dispersion in the extinction curves. It 
is evident that the MEC deviates substantially from the SMC, 
displaying a flattening at short wavelengths (Ar F < 2000 A), 
although the dispersion is high. 

Secondly, we perform the simultaneous fit of all the seven 
reddened quasars; we report the results of this second approach 
in Table [3] The resulting global best-fit extinction curve (GEC) 
is indicated by a black dashed line in Fig. |5J and is very close 
to the MEC. In this case, we can also infer the confidence lim- 
its, which are similar to those obtained for individual quasars, as 
shown by the hatched area in Fig. [5] Although the simultaneous 
fitting should improve and decrease the width of the confidence 
limits (thanks to the higher quality statistics than for individ- 
ual spectra), the confidence interval remains large implying that 
the curve of the SMC is consistent with the GEC. As we show 
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Fig. 2. Observed/rebinned spectra are indicated by black lines/filled blue squares, while red dashed lines represent the best fit spectra 
obtained through the^- 2 analysis. The gray bands report regions of poor atmospheric transmission. 



in the following, most of the dispersion and the width of the 
confidence intervals is due to the non-BAL quasars. Originally, 
the difference between BAL and non-BAL quasars was mostly 
ascribed to dif ferent viewing angles of objects that are intrin- 
sically similar ( Wevman n et alj|199lt iGoodrich & Millerl 19951 
iHines & Wills lll995b . However, the finding of systemic differ- 
ences between BAL and non-BAL quasars in terms of obser- 
vational properties that should be isotropic, has suggested that 
these two classes of quasars are intrinsically different. In par- 
ticular, the finding that the radio morphology of BAL is gen- 
erally more compact than those of non-BAL qua sars suggests 
that th e former are in an earlier evolutionary stage dBecker et al.l 
2000). If BAL and non-BAL quasars were to trace different 
evolutionary stages, this could also imply different dust prop- 
erties. Different extinction curve properties are already apparent 
in Fig. [4] the extinction curves of BAL quasars (solid curves) 
are systematically flatter than those of non-BAL quasars (dashed 
curves). In Fig. [6] we show the MEC and GEC for BAL and non- 
BAL quasars separately. The average extinction curve (MEC) of 
BAL quasars exhibits a much lower dispersion than the whole 
sample of reddened quasar. Moreover, for this subsample the 
confidence intervals on the simultaneous/global fit (GEC) are 
much narrower. As a consequence, for BAL quasars, the devi- 



ation of the simultaneous/global extinction curve from the SMC 
is far more significant than in the case of the whole sample of 
quasars; in particular, the deviation from the SMC is significant 
at a confidence level higher than 95%. The MEC and GEC of 
non-BAL quasars are intermediate between the extinction curve 
of BAL quasars and the SMC. The confidence intervals of the 
GEC for non-BAL quasars are much wider, making the GEC 
consistent with the SMC curve. 

Our speculative interpretation is that the extinction curve of 
high-z BAL quasars reflects the dust properties in the early evo- 
lutionary stages of their host galaxies, when SNe possibly pro- 
vided a preferential mechanism for dust production. High-z non- 
BAL quasars are probably in a more evolved stage, character- 
ized by dust properties that are in-between those of the SMC and 
those of primordial dust grains. However, we cannot exclude that 
the different extinction curve of high-z BAL quasars is related to 
a different reprocessing of dust grains in the ISM, from those of 
lower redshift and non-BAL quasars. For instance, higher den- 
sities may favor the growth and coagulation into larger grains 
flattening the extinction curve. The preferential destruction of 
small grains due to sublimation is also important; indeed the 
sublim ation radius for small grains is larger than that of large 
grains dLaor & Draine 1119931) . At high quasar luminosities, such 
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Fig. 3. Same as Fig. [2] We note a prominent Fell bump in the ULAS J0203+0012 spectrum. 



as those characterizing high-z quasars (because of selection ef- 
fect) this differential sublimation effect may play a major role. 
If the latter is the main explanation of our findings, then the 
extinction curve depends mainly not on redshift but on quasar 
luminosities. In this case, the flattening is only associated with 
high-z quasars and not with early galaxies. 

The higher confidence with which the GEC of BAL quasar 
deviates from the SMC is also most likely caused by observa- 
tional effects and in particular by BAL quasars being generally 
more reddened, and their extinction curve therefore being able 
to be measured with higher accuracy. 

We finally emphasize that our result is based on the analysis 
of individual quasars and not on the use of composite spectra, a 
method that can inf er biased (artificially flatter) extinction laws 
(IWillott et al.ll2005l) . 

3.5. Selection effects and spurious correlations 

Fig. [7] (left-hand panel) shows the extinction A3000 as a func- 
tion of the spectral index ax- The extinction of quasars that 
are not reddened according to our criteria is set to be zero. 
The empty triangle shows the average slope of "unreddened" 
quasars, (ax) = -1.68, which is fully consistent with the aver- 
age slope of unreddened quasars at low-intermediate redshifts 



(Reichard et al. 2003). This indicates that there is little or no 
evolution in the intrinsic spectrum of quasars at high-z and also 
supports our assumption that low-z templates are appropriate for 
modeling high-z quasars. 

However, we note that all quasars identified as "reddened" 
are characterized by very blue slopes, on the tail of low- 
intermediate quasar distribution. This is most likely the com- 
bined consequence of selection effects and our (conservative) 
criteria for identifying reddened quasars, as discussed in the fol- 
lowing. A reddened power law can still be fitted with a (redder) 
power law, and classified as "unreddened", as long as the extinc- 
tion is not too high. The deviations from a power law become 
significant only beyond a "minimum" extinction. For each in- 
trinsic ax, we inferred the minimum A3ooo,min for which the de- 
viations from a (redder) power law become significant and cause 
the target to be classified as "dust reddened". We plot A3ooo,min 
in Fig.|7]with a thick, black solid line, which has a value of about 
0.8 mag and is independent of the intrinsic a\. Quasars affected 
by extinction lower than this minimum value (gray shaded area) 
are not classified as "absorbed" according to our criteria. That 
A3ooo,min decreases at ax > -1.2 is simply a consequence of our 
limited range of intrinsic slopes ([-2.9; -0.2], Sect. I3.lt . which 
makes an absorbed quasar spectrum intrinsically red and difficult 
to fit with a redder slope, simply because of the "lack" of redder 
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Fig. 4. Best fit extinction curves of reddened quasars. The solid lines are for BAL quasars, while dashed lines are for non-BAL 
quasars. For comparison, the SMC extinction curve is also shown and labeled in the figure (dotted black line). The panel on the left 
shows the results assuming a minimum intrinsic slope a^ m i n = -2.9, while the panel on the right is obtained with ax, m in = -2.6. 
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Fig. 5. The solid line shows the mean extinction curve (MEC) computed by averaging the results obtained for all reddened quasars, 
while the shaded region shows the dispersion. The dashed line shows the extinction curve obtained by the simultaneous/global fit of 
all the reddened quasars (GEC), while the hatched area shows the 68% confidence limit. 
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slopes in the template library. If we remove the lower limit to 
the allowed intrinsic slopes, then we obtain the horizontal solid 
thick line, which extends the A3ooo,min ~ 0.8 line to ax = -0.3. 

Another selection effect is related to the parent sample of 
high-z quasars, which were selected as color dropouts (i-z at 
z~6 and r-i at z~4.5, which sample the IGM Lya forest break 
at the respective redshift), and to also have a blue color on the 
filters rewards of the dropouts (z-J and i-z). The latter selection 
criterion is required in order to minimize the contamination by 
brown dwarfs and low redshift interlopers. The latter require- 
ment can be translated into a maximum observed spectral slope 
(ci,max)- This implies a maximum extinction A3ooo,max affecting 
quasars, which corresponds to the reddening that causes the ob- 
served slope to be redder than a a, max- Clearly A3ooo,max depends 
on the intrinsic slope ax of the quasar, and it is shown with a 
dashed line in Fig. [71 where we have assumed au.max = -0.2, as 
inferred by the minimum of the slope distributions in the SDSS 
quasar survey. Quasars reddened by an extinction higher than 



the dashed line (gray hatched region) are lost by the quasar color 
selection criteria. 

Quasars surviving both our definition of reddened quasars 
(above the gray shaded area) and the survey color-selection cri- 
teria (below the gray hatched area) are clearly limited to a small 
region of the A3000 - cex domain, mostly in the low ax range, 
which is indeed where our reddened quasars are identified. 

Another selection effect is linked to the detectability of 
quasars as a function of spectral slope and extinction. For a given 
bolometric luminosity, quasars with a bluer slope are easier to 
detect because a higher fraction of the flux is emitted in the UV, 
which is the band where quasars are detected. As a consequence, 
for a given bolometric luminosity, dust extinction brings intrin- 
sically red quasars below the survey detection threshold more 
easily than blue quasars. This effect is illustrated by the dotted 
lines in Fig. [71 which show the maximum A3000 at which a quasar 
with a given bolometric luminosity can still be detected by the 
limiting magnitude of the survey in the i-band (z ~ 4, green 
lines) or in the z-band (z ~ 6, orange lines) as a function of the 
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Fig. 6. Mean (MEC) and simultaneous/global (GEC) extinction curve of reddened quasars divided into BAL (left) and non-BAL 
(right). The coding is the same as in Fig. [5] 
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Fig. 7. Left hand panel: Red asterisks show the best fitting values of A3000 and ax, as obtained for reddened quasars. Red asterisks 
surrounded by blue squares represent BAL quasars. The black filled square represents the corresponding mean value, and the error 
bars the standard deviations. The empty black triangle and the associated error bars denote the ax mean value and the lcr dispersion 
in the case of unreddened quasars. Solid orange and dotted green curves represent the maximum value of A3000, such that, for a 
given bolometric luminosity Lb„i, the observed flux at a given wavelength is greater than the flux detection limit in the Z-band and 
I-band, respectively. The bolometric luminosity has been varied in the range 10 47 < Lboi < 10 48 (curves from the bottom to the top). 
The orizontal black solid line represents the minimum extinction required for a quasar to be classified as "dust reddened". Quasar 
spectra reddened by an extinction lower than this minimum value (gray shaded area) can still be fitted with a power law redder 
than the intrinsic one, and classified as "unreddened". The black dashed line shows the effect of quasar color selection criteria on 
our results. High-z quasars are selected to have a blue color on the filters redwards of the dropout (z-J and i-z). This requirement 
can be translated into a maximum observed spectral slope (a^max) which implies a maximum extinction A3ooo,max affecting quasars. 
Clearly A3ooo,max depends on the intrinsic slope ax of the quasar, and it is shown with the black dashed line for au jmax = -0.2, i.e. the 
minimum of the slope distributions in the SDSS quasar survey. Quasars reddened by an extinction higher than the dashed line (gray 
hatched region) are lost by the quasar color selection criteria. Right hand panel: Trend of A3000 with redshift. The lines represent the 
maximum value of A3000 for which a reddened spectrum, for a fixed bolometric luminosity (Li, i = 10 48 erg s" 1 ), is still detectable 
by the SDSS. The results have been obtained assuming a slope for the intrinsic continuum ax = -2.4 (dotted), ax = -1.6 (dashed). 
The horizontal black solid line and the gray shaded region have the same meaning of the ones in the left hand panel. 



spectrum slope. For each of the two cases, the three lines corre- 
spond to different luminosities (log Lboi = 47,47.5,48 erg s~'). 
This selection effect clearly also contributes to bias the sample of 
reddened quasars against intrinsically red objects. However, the 
true magnitude of this effect is difficult to quantify because of 
claims of a possible correlation between luminosity and spectral 



slope in quasars (e.g. Pu et al. 2006). As a consequence, the de- 
pendence of the maximum A3000 on ax may be globally steeper 
for the real population of high-z quasars than illustrated by the 
dotted lines in the left hand panel of Fig. [7] 

Figure [7] (right hand panel) shows the extinction A3000 as a 
function of redshift for reddened quasars, showing an anticor- 
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relation between these two quantities. However, this anticorre- 
lation is also probably mostly driven by selection effects: for a 
given bolometric luminosity, more distant quasars appear fainter 
and therefore fall below the detection threshold of the survey 
with a lower amount of dust extinction. To illustrate this effect, 
the dashed and dotted lines in Fig.|7](right hand panel) illustrate 
the maximum extinction allowed for quasars with bolometric lu- 
minosity of 10 48 ergs s _1 (with a = -2.4 and a = -1.6, re- 
spectively) in order not to be missed by the SDSS survey, which 
approximately reproduces the trend observed in our reddened 
quasars. The "jump" at z~5 is caused by the different photo- 
metric bands used to select quasars in the two redshift ranges. 
Our reddened quasars approximately follow the expected rela- 
tion. Although our data have very limited statistics, at very high 
redshift (z>6.3) it seems that there is a shortage of highly red- 
dened quasars that are still expected within the selection effects 
(extrapolation of the dashed and dotted curves). The latter re- 
sult may be indicative of a real lack of dust at these very early 
epochs because of the shortage of ti me to produce dust, in line 
with the finding of Jian g et alj d2010h . Another possible effect is 
that, since the Lya cutoff of these objects is partly outside the 
z-band, they would need to have very high instrinsic luminosi- 
ties to be reddened and observed by SDSS; since the luminosity 
function is very steep at the bright end, this may cause a short- 
age of reddened quasars at z> 6.3 relative to their lower redshift 
counterparts. However, the statistical significance of our result 
is still marginal. Finally, we note that there is a lack of quasars 
at z<5 with extinctions A3000 in the range 0.8-1.8 mag, which 
should be in principle observed. The statistics are low enough to 
ensure that this shortage of quasars with intermediate extinction 
is not a real worry. However, we note that at z<5 many quasars 
are only observed with FORS2, hence their data have a very 
limited spectral coverage; as a consequence, for these quasars 
the minimum extinction required to be classified as "reddened" 
is significantly higher, approaching A3ooo,min ~ 1.5 mag, hence 
partially explaining the additional selection effect noted in Fig. [7] 

3.6. Degeneracy between extinction curves 

An SMC-extincted spectrum with a given slope is very similar to 
an intrinsically redder spectrum absorbed by an extinction curve 
that, relative to the SMC, flattens in the UV. This effect is illus- 
trated in Fig. [8] which shows two spectra with intrinsic slopes 
of a,\ — -1.6 and ax = -2.0, but absorbed with the MEC and 
SMC extinction curves, respectively. The two resulting spectra 
are indeed similar, although we verified by means of simulations 
that with our analysis we can recover the correct extinction curve 
in both cases. However, what is interesting within this context, 
is that any bias introduced by this degeneracy goes in the op- 
posite direction of the results obtained by us. Indeed, that we 
find a deviation from the SMC extinction curve flattening in the 
UV (the MEC), and a blue intrinsic slope (as discussed in the 
previous section) strengthens our finding, meaning that in an un- 
biased quasar sample (i.e. with a lower limiting flux, hence with 
an intrinsic slope closer to the average of unreddened quasars) 
the deviation from the SMC may be even more pronounced than 
the conservative result reported here. 

We finally note that there is an unavoidable degeneracy in 
the slope of the extinction curve, which is intrinsic in any study 
using only the spectral shape (i.e. without knowledge about the 
intrinsic luminosity of the object). Given an extinction curve 
tj = A,i/A3ooo> fitting an observed spectral shape, any other ex- 
tinction curve of the form fC = k(f^ - 1) + 1, can indeed fit 
equally well the same observed spectrum with any value of k, by 
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simply changing the absolute extinction A3ooo- We reiterate that 
this degeneracy affects any study attempting to infer extinction 
curves solely based on the objects spectral shape (i.e. the vast 
majority of previous studies). Only in cases where the intrinsic 
luminosity is known, and therefore A3000 is not free to vary, can 
the slope (k) of the extinction law be determined. The alterna- 
tive approach, as adopted here, is to trust only a limited library 
of (physically justified) extinction curves, without leaving the 
possibility of changing their individual slope. 

4. Implications for primordial galaxies 

4.1. Observational implications 

We have inferred the extinction curve for quasars at z>4. We pro- 
pose that the same extinction curve applies also to star-forming 
galaxies at high-z in general and adopt this as a working assump- 
tion. If the dust properties at high redshift differed significantly 
from those at low redshift, this would have important effects for 
the observational properties of high-z galaxies. The extinction 
curve adopted to correct the rest-frame UV radiation of high-z 
galaxies should be modified with respect to low redshift tem- 
plates, otherwise the inferred intrinsic luminosity and spectral 
shape (used to infer the star-formation rate and stellar popula- 
tion) will be erroneously estimated. In particular, the evolution 
of the cosmic star-formation rate at z > 4 may be erroneously 
traced if the change in the dust extinction properties is not prop- 
erly accounted for. In this section, we discuss, and when possible 
quantify, these observational implications. 

4.1 .1 . Attenuation curves 

For quasars, the simple dusty "screen" geometry applies, but for 
galaxies one has to consider that dust is mixed with the emit- 
ting sources (either stars or ionized gas). As a consequence, for 
galaxies it is more practical to use "attenuation curves", i.e. the 
ratio of the observed spectrum to the intrinsic total light emitted 
by the whole system, before dust absorption. The attenuation 
curve is of course associated with the extinction curve (which is 
intrinsic to the dust properties), but also depends on the geome- 
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Fig. 9. Left hand panel: The attenuation curve corresponding to the SMC extinction curve (solid magenta line) is shown by an 
orange dashed line. For comparison, the Calzetti law is also indicated by a dotted green line. It has been shown that the Calzetti law 
is the result of the attenuation by a medium with an SMC-type dust (IGordon et al.lll997t llnoue II2005I) . We can recover this result 
using Eq. 3. Right hand panel: Same as in the left hand panel, but for the MEC. The attenuation curve resulting from the MEC is 
flatter than the Calzetti law, which is generally adopted for quantifying the reddening of galaxies at all epochs. 



try of the emitting sources relative to the dusty medium, as well 
as the viewing angle. We attempt to derive a simple version of 
the attenuation curve associated with the MEC extinction curve, 
which can be applied to high-z galaxies. However, we also con- 
sider the simple "screen" cas e, since this geom etry seems to ap- 
ply to some high-z galaxies (Red dv et al.ll2.Q10b - 

For galaxies in the local Universe, the effect of dust is as- 
cribed to the attenuation law by Calzetti et al. (1994). It has been 
shown that the Calzetti law is the result of the attenua t ion by 
a med ium with an SMC-type dust (IGordon et al.lll997t llnoue I 
120051) . Therefore, that the MEC is flatter in the UV relative to the 
SMC indicates that the attenuation law resulting from a MEC- 
type dust must be flatter than the Calzetti one. We assume a very 
simple case where dust and stars are uniformly mixed and the 
density is constant throughout the galaxy. In this case, it is easy 
to infer that the attenuation curve normalized to 3000 A is given 
by 

A a„ l °8io \^r- 

a a _ \ Sio 73000 / (3) 
A 3ooo k*io(-S5rJ 

where A^/Aaooo is the dust extinction curve, and T3000 is the total 
dust optical depth at 3000 A of the dust throughout the galaxy 
(i.e. the dust optical depth integrated over the entire thickness of 
the galaxy along the line of sight). First, we insert in Eq. (O the 
SMC extinction curve, and we show in Fig.[9]that we can almost 
perfectly reproduce the Calzetti law with T3000 = 6. We then in- 
vestigate the effect of changing the underlying extinction curve 
by keeping the same value of T3000 an d using the MEC instead 
of the SMC; we refer to the resulting attenuation curve with the 
notation MEC att . The result is shown in the right-hand panel of 
Fig. [9] where the blue solid line denotes the MEC, the dashed red 
line the MEC fl », and the green dotted line the Calzetti law, which 
is reported for comparison. We note that for this analysis we use 
the MEC obtained from the analysis of BAL quasars, which (as 
discussed in Sect. 13. 4t we consider as the truly "primordial" ex- 
tinction curve. We indeed confirm our expectations, i.e. that the 



attenuation curve inferred by the MEC is flatter than the Calzetti 
one. 

4.1 .2. Reddened galaxy spectra 

We now investigate how different extinction and attenuation 
curves affect the shape of a typical galaxy spectrum. As an exam- 
ple, we consider a (unabsorb ed) spectrum of a STARBURST99 
model (Leithe rer et alJ[l999l) for a galaxy 70 Myr old that has 
experienced an instantaneous burst of star-formation, resulting 
in the production of 10 6 M Q of stars. We assume a metallicity of 
Z = 5 x 10~ 2 Z Q and apply the constraint that the spectral region 
blueward of the Lya is completely absorbed by the intervening 
intergalactic medium (which is appropriate for galaxies at red- 
shifts approaching the reionization epoch). Since we only aim to 
make an illustrative comparison between the extinction effects of 
SMC and MEC, we neglect the effect of nebular lines, although 
they are certainly promine nt in early star-forming gala xies, es- 
pecially at low metallicity dSchaerer & de Barros II2010I) . 

For the case of pure screen extinction, in the left-hand panel 
of Fig. [10] we compare the effect of the SMC-like (magenta 
dotted line) and MEC-like (blue dashed line) reddening on the 
spectrum of the high-z star-forming galaxy template (black solid 
line), by assuming the same extinction at 3000 A (A3000 = 
2 mag). 

The MEC clearly makes the spectrum bluer and more lumi- 
nous than the SMC, as expected. This first result is particularly 
propitious for the knowledge of the high-z Universe, since it im- 
plies a much easier detection of early galaxies with future facili- 
ties (e.g. JWST, ELT). In the specific case considered here, there 
is a difference of nearly a factor of ten between continuum inten- 
sity at the Ly-break between the use of the two different extinc- 
tion curves. What is most relevant is that if the apparently blue 
spectrum produced by the MEC extinction is interpreted using 
the SMC extinction curve, then one infers that the dust extinc- 
tion is much lower than the one actually affecting the spectrum. 

In the right-hand panel of Fig.QJJ we show the reddening of 
galaxy spectra using the attenuation curves and, in particular, the 
Calzetti law versus the MEC attenuation curve. In this case, the 
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Fig. 10. Spectra of galaxies at redshifts higher than the reionization epoch. The solid black lines represent unabsorbed templates for 
a galaxy 70 Myr old, having a metallicity Z = 5 x 1O~ 2 Z , which has experienced an instantaneous star-formation, resulting in the 
production of 10 6 M of stars. Dotted and dashed lines show the effects of different types of reddening, as labeled in the figure. The 
spectral region blueward of the Lya is assumed to be completely absorbed by the intervening intergalactic medium. 



effect is much lower than in the case of pure screen extinction. 
However, in the case of attenuation by MEC the spectrum still 
appears bluer than in the case of attenuation with the Calzetti 
law, hence leading to an underestimate of the effect of dust if not 
interpreted properly. We quantify this effect in the next section. 

4.1 .3. Correction for dust extinction in galaxies 

For high-z galaxies, which are generally much fainter than 
quasars, only a few broad-band photometric measurements of 
the rest-frame UV continuum are usually available. As a conse- 
quence, the effect of dust reddening is often computed by simply 
deriving the spectral slope /? (F cont oc A 13 ), and measuring its devi- 
ation from an average slope assumed to represent the unabsorbed 
case (/? ~ -2.2). The relationship between slope /? and extinc- 
tion at Atest = 1600 A was calibrated by iMeurer et al.l d 19991) 
and subsequently adapted by iBouwens et alf i 20091) tothe case 
of galaxies at high-z, whose slopes are inferre d by using only 
two b road-band filters. The relation inferred by IBouwens et al.l 
(120091) is 



A 1600 = 4.91 +2.21/?, 



(4) 



where /? is computed in the wavelength range between 1600 A 
and 2300 A. The resulting relation is shown in Fig. [TT] with a 
black solid line. The relation between extinction and spectral 
slope obviously depends on the extinction or attenuation curve 
assumed. The relation in Eq. [4] can t> e reproduced in terms of 
extinction by a Calzetti attenuation curve, as illustrated by the 
green dotted line in Fig. QTJ (see also Table 2 of Meurer et al. 
1999). When we adopt the attenuation curve associated with the 
MEC at z>4, then the relation between A1600 and /? is shallower, 
as a consequence of the attenuation curve being flatter in the UV 
than the Calzetti curve. More specifically, in this case, the rela- 
tion given by Eq.|4]becomes 



A 



1600 



9.61+4.33/?, 



(5) 



where we keep the same assumption about the average intrinsic 
slope as in Eq.|U i.e./? ~ -2.2. We recall that Eq.|5]is obtained by 
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Fig. 11. Slope (/?) of the UV-continuum of a galaxy as a function 
of the extinction/attenuation at 1600 A (Ai6oo)- The solid black 
line show the Meurer et al. 1999 relation. The green dotted (red 
dashed) lines shows the /?-Ai6oo relation for a galaxy with an 
intrinsic slope of /? = -2.2, assuming that the dust attenuation is 
described by the Calzetti law (MEC fl „). 



adopting the MEC attenuation curve, computed based on the as- 
sumption that dust and stars are uniformly mixed and the density 
is constant through th e galaxy. The latter re lation is compared 
with the one used by IBouwens et al.1 (120091) in Fig. [TT] If the 
reddening curve inferred for quasars at z>4 also applies to high- 
z galaxies, then we suggest that Eq. [5] is more appropriate for 
describing the extinction in high-z galaxies than Eq. [4] Both of 
them of course rely on the assumption that the average intrinsic 
slope /? is about -2. 2, which may not be appropr iate for high-z 
metal poor galaxies (Schaere r & de Barros I i2010). although our 
main goal here is simply to show the effect of the different ex- 
tinction curve on the inferred correction of the luminosity (and 
star-formation rate) in high-z galaxies. 
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Fig. 12. Effects of different dust corrections on the star-formation rate evolution for luminous (on the left) and less luminous (on 
the right) galaxies. Black asterisks and filled green squares represent, respectively, the observed and corrected star-formation rate 
(along with the corresponding errors) by Bouwens et al. (2009), assuming the ft - Ai^oo relation by Meurer et al. 1999. The red 
dashed lines report the corrected star-formation rate, obtained taking into account the results of this study. 



4.1.4. The star-formation rate at high-z 

Dust correction is a key requirement when measuring the star- 
formation rate density evolution. Observations of high-z galax- 
ies are consistent with a star-formation rate density at z ~ 6 that 
is ~ 6 times less than that inferred for lyman break galaxies at 
-3-4 (Bunker et al. 2004; Bunker et al. 2009; Bouwens et 
al. 2009). At z~7, the star-formation rate density seems to be 
even lower, a factor of about 10 less than the one inferred for 
low-z LBGs. Moreover. [Bouwens etal] d2009l) find that the UV- 
continuum slope of z ~ 5 - 6 galaxies is substantially bluer than 
Z ~ 2-4 ones, therefore claiming that high-z galaxies are less ob- 
scured by dust than their lower redshift counterparts. However, it 
is intriguing that the change in slope of the continuum of galax- 
ies occurs just at the redshift (z~4-5) where we find that the ex- 
tinction curve may change. This suggests that the change in the 
average continuum slope of galaxies is actually due to a change 
in the dust properties responsible for reddening the spectra. 

iBouwens et alj d2009f) indeed mention that changes in the 
mean dust extinction curve may be responsible for the redshift 
dependence of the UV-continuum slope, although their dust cor- 
rections, used to infer the star-formation rate density at high- 
z, are computed starting from empirical relation between A1600 
and /3 obtained for galaxies in the local Universe (Eq.|U). In the 
following, we investigate the effect on the star-formation rate 
density when the redshift dependence of the average continuum 
slope is ascribed to a variation in the dust extinction curve at 
high-z. 

We do not aim to perform a very detailed investigation, 
which would require a remodelling of all the photometric data 
points in galaxies at z>4 with the MEC attenuation curve, which 
will instead be the topic of a separate, dedicated paper. We sim- 
pl y take the star-formatio n rate density observationally inferred 
bv lBouwens et alj (|2009) and recompute the dust extinction cor- 
rections by using the relation between A1600 and /3 given by 
Eq.|5] which is appropriate for the MEC attenuation curve. We 
assume that the star-formation rate is simply proportional to the 
UV luminosity at 1600 A. Figure fT2l shows the star-fo rmation 
rate density at z>2 inferred by Bou wens et alj d2009l) uncor- 
rected for dust extinction (black asterisks) and corrected accord- 



ing to Eq.[4] (green squares), i.e. by assuming the local relation 
between extinction and spectral slope. The left hand panel is for 
luminous sources (L > 0.3 L* =3 ), while the right hand panel is 
for less luminous sources (L > 0.04 L* =3 ). The red hatched area 
shows the star-formation rate corrected by using Eq. [5] (i.e. by 
assuming the MEC attenuation curve) and the observed spec- 
tral slopes given by IBouwens et alj (|2009). The width of the 
hatched area reflects the uncertainties in the spectral slopes. As 
expected, the star-formation rate density after applying the cor- 
rection obtained by using the MEC attenuation curve is substan- 
tially higher than previous estimates based on local calibrations 
of the dust absorption correction. With the MEC-attenuation cor- 
rection, the decline in the SFR density is less pronounced than 
previous estimates. 

We note that if the screen appr oximation holds fo r young 
high-z sources (as suggested by iReddv et alJ 1201 Oh . which 
are possibly associated with dusty winds acting as foreground 
screens, then the magnitude of the differential correction on the 
SFR density would be much larger than estimated above. In the 
case of foreground screen-like absorption, one should indeed use 
the extinction curves (and not the attenuation curves), and the 
difference between MEC and SMC should be much larger than 
for the corresponding attenuation curves. 

We caution that we do not intend to provide a new measure- 
ment of the SFR at z > 4, but simply to show how the actual SFR 
may be underestimated by adopting a relationship between spec- 
tral slope and extinction inferred for local or intermediate red- 
shift galaxies (Eq.[4]i. As mentioned above, the careful derivation 
of the SFR by using a different prescription for the dust extinc- 
tion requires a new, dedicated modeling of individual galaxies at 
z>4, possibly by including additional photometric (or even spec- 
troscopic) measurements, which will probably become possible 
only with the next generation of facilities (JWST, TMT, EELT). 
Here, we note that our results may skew the peak of the SFR den- 
sity from z ~ 2 - 3 (Hopkins & Beacom 2006) towards higher 
redshifts. 

Our results may imply that the sources responsible for the 
reionization (or for keeping the universe reionized), and cur- 
rently not detected in the high-z surveys, are not below the de- 
tection limit because they are intrinsically faint, but because they 
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have been absorbed by dust. The knowledge of dust properties in 
the early Universe is therefore funda mental to studies of high-z 
galaxies during reionization (see also lDaval et al .1120 1 Oalfbb . 

5. Conclusions 

We have investigated the extinction properties of dust at 
3.9<z<6.4 by using optical-near infrared spectra of 33 quasars 
in this redshift range. Our main goal has been to investigate 
whether the SMC extinction curve, which accurately describes 
quasar reddening at z < 2.2, also represents a good prescription 
at these epochs. 

We have fitted the quasar spectra with templates by varying 
the intrinsic spectral slope and by using a large grid of extinction 
curves, including the well-known empirical laws and several the- 
oretical extinction curves (as well as mixtures of these extinction 
curves with the SMC). We find that 7 quasars in our sample re- 
quire substantial (0.8 > A3000 ^ 2.0) extinction. The best-fit 
extinction curve was found to consistently differ from the SMC, 
although the latter is still (marginally) consistent with the indi- 
vidual spectra. The best-fit extinction curve is generally flatter 
than the SMC in the UV, although there is a large spread among 
the different extinction curves. We note that most of the disper- 
sion among the best-fit extinction curves originate in the mixture 
of BAL and non-BAL quasars. BAL quasars are always charac- 
terized by flatter extinction curves relative to non-BAL quasars. 
If BAL reddened quasars are fitted simultaneously then the de- 
viation from the SMC becomes highly significant (> 2<x). The 
extinction curves characterizing non-BAL quasars have instead 
properties that are intermediate between BAL quasars and the 
SMC, and the deviation from the latter is much less significant. 
We suggest that reddened BAL quasars at high-z are in an earlier 
evolutionary stage (in agreement with previous studies) when the 
first generation of dust is produced; then, as quasars evolve into 
non-BAL reddened systems their dust extinction curve proper- 
ties approach the SMC extinction. However, the higher confi- 
dence with which BAL quasars deviate from the SMC is proba- 
bly also caused by BAL quasars in general being more reddened 
and therefore their extinction curve being able to be measured 
with higher accuracy. 

We infer a mean extinction curve (MEC) for all reddened 
quasars at z>4, and by also dividing them into BAL and non- 
BAL. According to the scenario discussed above, the MEC of 
BAL quasars should represent the most appropriate description 
of the extinction curve in early Universe galaxies. 

The finding that the extinction curves at z>4 deviate from 
those observed at lower redshift suggests that the dominant dust 
production mechanism at these early epochs differs from that 
of local or intermediate redshift objects. At these early epochs, 
there is probably an higher contribution of dust by SNe than 
AGB stars. This is in line with the extinction curves theoreti- 
cally expected for SN-like dust providing generally the best fit 
to t he quasar spec t ra at z >4. 

IValiante et al. (2009) demonstrated that the relative contri- 
bution of AGBs and SNe to the dust production in quasar is 
strongly dependent on the star-formation history of the host 
galaxy. In particular, the SNe contribution dominantes until the 
last major burst of star-formation. This scenario is supported by 
our finding of a more pronounced deviation from the SMC in 
the case of BAL quasars. In fact, since BAL quasars are proba- 
bly younger than no-BAL quasars, the former may be closer to 
the starburst phase than the latter. 

We suggest that the extinction curve inferred for quasars at 
z>4 may also apply to young galaxies at similar redshifts. More 



specifically, we suggest that the average change in the spectral 
slope for galaxies at z>4 (getting bluer) may reflect the same 
variation in extinction curve observed in quasars, rather than a 
lower dust content. We infer the attenuation curve associated 
with the MEC of BAL quasars. The relationship between galaxy 
spectral slope f3 and extinction A1600 used for local and inter- 
mediate redshift galaxies is modified using the MEC attenuation 
curve, which we regard to be more appropriate at z>4. We use 
this relation to correct the star-formation rate density at z~4-6 
for dust extinction, leading to substantially higher estimates than 
derived previously. 
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Table A.l. Unreddened quasars 



Name 


z 


&A,Noext 




SDSS J1030 


6.28 


-2.05 


1.00 


CFHQS 1509 a 


6.12 


-1.29 


0.71 


SDSS J 1602 


6.07 


-0.82 


1.00 


SDSS J1630 


6.06 


-1.71 


0.08 


SDSS J 1306 


5.99 


-2.27 


0.19 


SDSS J1411 


5.93 


-1.66 


0.88 


SDSS J0836 


5.80 


-1.87 


0.97 


SDSS J0005 


5.85 


-1.42 


1.00 


SDSS J0002 


5.80 


-1.60 


0.06 


SDSS J1044" 


5.78 


-1.65 


0.32 


SDSS J02313 


5.41 


-2.71 


0.33 


SDSS J1614 


5.31 


-1.22 


1.00 


SDSSpJ1208 


5.27 


-2.62 


0.59 


SDSS J0957 


5.16 


-2.66 


0.36 


SDSS J0756 b 


5.05 


-1.64 


0.98 


SDSS J 1204 


5.05 


-2.10 


0.24 


SDSS J0017" 


5.01 


-1.95 


0.10 


SDSS J0338 


5.00 


-1.87 


0.42 


SDSS J2216 


4.99 


-1.49 


0.45 


SDSS J1605 ac 


4.92 


-1.33 


1.00 


SDSS J2200 


4.77 


-1.37 


0.48 


SDSS J0120 a 


4.73 


-1.21 


1.00 


SDSSp J1021 


4.70 


-1.01 


0.86 


SDSS J 1603 


4.39 


-1.78 


0.82 


SDSS J0156" 


4.32 


-1.22 


0.91 


SDSS J0239 a 


4.02 


-1.39 


1.00 



a BAL quasars; 

b For this spectrum we find a P(x~) ~ 0.6 in the case of reddening by 

an extinction curve given by 30% SMC + 70% PISN2; 
c For this spectrum we find a Pix 1 ) ~ 0.6 in the case of reddening by 
an extinction curve given by 70% SMC + 30% H2; 



tion curves is very extended, we plan to reobserve these object 
to improve our analysis. 



Appendix A: Results for unreddened quasars 

In Fig lA.H4A.5l we show the spectra of all the unreddened 
quasars (black solid line) along with their best fitting template 
(red dashed line). As for the reddened quasars, blue squares 
show the rebinned spectral value used for the fitting procedure. 
In Table I A. fl we report the fitting results for quasars that do not 
require substantial reddening. We note that there are some ob- 
jects for which the "no extinction" case does not provide a good 
fit, although they are considered unreddened, since the best-fit 
extinction curve is also characterized by a P(Xbf > 68%). This 
may be due to both the actual extinction curve not being part of 
our grid and some problems in the data. Since the grid of extinc- 
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A obs [>m] " ' X ob S [M m ] 

1.0 1.2 1.4 1.6 1.8 2.0 2.2 1.0 1.2 1.4 1.6 1.8 2.0 2.2 




1500 2000 2500 3000 1500 2000 2500 3000 

X RF [A] X RF [A] 



Fig. A.l. Unreddened quasar spectra. The observed/rebinned spectra are indicated by black lines/filled blue squares, while the red 
dashed lines represent the best fitting spectra resulting from the y 1 analysis. The gray bands report regions of bad atmospheric 
transmission. In each panel, the name and the redshift of the quasar is shown at the top left, while the best fitting slope is reported 
at the top right. 
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Fig. A.2. Fig. A.l continued. 
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Fig. A.3. Fig. A.l continued. 



1200 1300 1400 1500 1600 1700 1800 

V[A] 



19 



S. Gallerani et al.: The extinction law at high redshift and its implications 
A obs [>m] " ' ^cbs [M m ] 

1.0 1.5 2.0 1.0 1.5 2.0 




1200 1300 1400 1500 1600 1700 1800 1500 2000 2500 3000 3500 4000 4500 

X RF [A] X RF [A] 

. A.4. Fig. A.l continued. 
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